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Abstract 



The r decays to six-pion final states have been studied with the CLEO de- 
tector at the Cornell Electron Storage Ring. The measured branching frac- 
tions are B(t~ -» 2ir-ir+3ir°v T ) = (2.2 ± 0.3 ± 0.4) x 10" 4 and B(t~ -> 
= (1.7±0.2±0.2) x 10 4 . A search for substructure in these de- 
cays shows that they are saturated by intermediate states with r\ or u mesons. 
We present the first observation of the decay r~ — > 2ir~ir + ujv T and the branch- 
ing fraction is measured to be (1.2±0.2±0.1) x 10 -4 . The measured branching 
fractions are in good agreement with the isospin expectations but somewhat 
below the Conserved- Vector-Current predictions. 
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The decays of the r lepton provide a good test of the Standard Model predictions for the 
hadronic weak current. The six-pion branching fractions are related to the isovector part 
of the e + e~ annihilation cross section by the Conserved Vector Current (CVC) hypothesis. 
Isospin symmetry relates the relative branching fractions of the decays r~ — > 27r~7r + 37r°z/ T , 
r~ — > 37r~27r + 7r°z/ r , and r~ — > 7r~57r°z/ T . Therefore, the study of six-pion decays can be used 
to test the CVC hypothesis and isospin predictions. A better understanding of the resonance 
substructure in the decays is also of particular interest because of the potential application 
in suppressing the hadronic background in the measurement of the r neutrino mass. In this 
Letter, we present a study of the decays [[!) t~ — > 27r~7r + 37r z/ r and r _ — > 3tt~27t + 7t°v t . This 
includes measurements of the branching fractions and the search for resonance substructure. 
The latter results are used to identify the vector and axial-vector current contributions to 
the inclusive branching fractions and allow the proper comparison of the results with the 
CVC and isospin symmetry predictions. 

The data used in this analysis were collected with the CLEO detector |2|] at the Cornell 
Electron Storage Ring (CESR) at center-of-mass energy of 10.6 GeV. The sample corresponds 
to a total integrated luminosity 13.5 fb _1 and contains 12.3 x 10 6 t + t~ events 0. 

For the decay r~ — > 27r~7r + 37T°z/ r (37r~27r + 7r°z/ T ) we select events with four (six) charged 
tracks and zero net charge. The momentum of each track must be greater than 100 MeV 
and the polar angle 6 of each track with respect to the beam must satisfy | cos^l < 0.90. 
The track must be consistent with originating from the e + e~ interaction point. This vertex 
requirement also suppresses the background from r events with a K$ or photon conversion. 
The Ks background is further reduced by rejecting events containing a detached vertex with 
a 7r + 7r~ mass consistent with the nominal Ks mass. 

We define two exclusive sets of photon candidates in the barrel calorimeter (|cos#| < 
0.80): high quality (HQ) and low quality (LQ) photons. The selection criteria for HQ photons 
are designed to minimize the contamination of fake showers from hadronic interactions, while 
the acceptance of LQ photons ensures high event detection efficiency A HQ photon must 
have an energy E 1 > 120 MeV and a lateral profile of energy deposition consistent with 
that expected for a photon. The Fisher discriminant technique [|J is used to differentiate 
a real photon from a fake photon. The discriminant function is a linear combination of 
the energy of the shower and its distance to the intersection of the nearest charged track 
with the calorimeter surface. We use Monte Carlo simulated events (see below) to obtain 
the discriminant function that maximizes the separation between real and fake photons. A 
LQ photon is defined as a shower that does not pass the HQ photon requirements but has 
Ej > 30 MeV. The Fisher discriminant for differentiating a real photon from a fake photon 
is only used in the decay t~ — > 27r~7r + 37r°z/ T due to the higher photon multiplicity. HQ 
photons in the endcap calorimeter (0.80 < | cos^l < 0.95) are selected as in the barrel. 

Each event is divided into two hemispheres using the plane perpendicular to the thrust 
axis |J , calculated using both charged tracks and photons. There must be one charged track 
in one hemisphere (tag) recoiling against three or five charged tracks in the other hemisphere 
(signal), depending on the decay mode. 

In the tag hemisphere, the total invariant mass of charged tracks and photons must 
satisfy M tag < 0.5 GeV. In the signal hemisphere, there should be at least six (two) photons 
forming three (one) ir° candidates for the decay r _ — > 27r~7r + 37r°z/ r (37r~27r + 7r°z/ r ). The 7r° 
candidates are reconstructed from photon pairs in the barrel calorimeter. All HQ photons in 
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the barrel must be used in the reconstruction and no HQ photon in the endcap is allowed. 
In case of multiple entries, we select the photon combination with smallest \ 2 — Y^i=x[S~Ji, 
where S" 77 = (m 77 — m^/a^ (<r 77 is the mass resolution calculated from the energy 
and angular resolution of each photon). The signal region for the ir° candidates is defined 
as —3.5 < Sjry < 2.5. In the case of r _ — > 3ti^2-k + ti q u t we use sideband subtraction 
to estimate the fake n° background, with sidebands defined as —8.5 < >S 77 < —5.5 and 
4.5 < S 7J < 7.5. The total invariant mass of the hadronic system in the signal hemisphere 
must satisfy M 6n < M T = 1.777 GeV 0. The signal hemisphere must have a positive 
pseudo-neutrino mass: 

M 2 M: ■ Ml -2M ; i: (i:: - , 

where E§ n is the energy of the six-pion system in the r rest frame, assuming that the r has 
the full beam energy and that the r direction is given by the momentum vector of the six- 
pion system. This cut selects events with tau-like kinematics, suppressing both the hadronic 
background and the r migration background from lower multiplicity decays where the six- 
pion momentum is not as good of an approximation of the r direction. The background from 
r decays with Dalitz tt° decays or photon conversion in inner detector material is suppressed 
by requiring that any pairs of oppositely charged tracks have invariant mass exceeding 120 
MeV when at least one of the pair is identified as an electron and assigning the electron 
mass to both tracks. The electron identification requires that a track is matched to a shower 
with shower energy to momentum ratio exceeding 0.85, and, if available, specific ionization 
in the drift chamber consistent with that expected for an electron. Two-photon backgrounds 
are eliminated by requiring the direction of the missing momentum of the event to satisfy 

| COS ^missing | ^ 0.9. 

The detection efficiency and r migration background are calculated with a Monte Carlo 
technique. We use the KORALB/TAUOLA program || for the r event simulation. The 
decay t~ — > 27r~7r + 37r°z/ r is modeled using a mixture of r~ — > 27T~7T + i]b> T , t~ — > -n~2^r\v T1 
and r~ — > -k^2ti q ujv t . The other decay t~ — > 37r~27r + 7r z/ r is modeled using a mixture 
of r~ — > 2n~'n- + r]v T and r~ — > 27r~7r + Ci;z/ T . The relative mixtures are determined from 
the measured branching fractions presented in this and a previous Letter |?[] . We assume 
that the 3irr] decays proceed through nfi with a spectral function dominated by the form 
factor of the ai(1260) resonance ||. The 3ttuj system is modeled assuming dominance of the 
p(1700) resonance. The detector response is simulated using the GEANT program 0. The 
hadronic background is calculated empirically using a sample of high-mass tagged events 
(1.8 < M tag < 2.8 GeV) obtained from the data with the assumption that M tag and M 6tt are 
not correlated. The invariant mass spectra of the six-pion hadronic systems after applying 
all selection criteria are shown in Fig. [I]. There is good agreement between the data and 
expectation. The signal, background, and detection efficiency are summarized in Table Q. 
The t background includes the decays with K$ — > 7r + 7r~ and 7r°7r°. 

The six-pion r decays can proceed through the rj or u intermediate hadronic states. 
Figure |2| shows the invariant mass spectra of 3ir combinations in the r] mass region. To 
improve the resolution, the ir° candidates have been kinematically constrained to the nominal 
7T° mass. In the search for 77 — > 7r + 7r~7r°, we reduce the combinatoric background by rejecting 
events with 3tt° invariant mass within 20 MeV (~3cr) of the nominal 77 mass. There are 
clear enhancements at the rj mass in all three spectra. The invariant mass spectra of 3n 
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TABLE I. Summary of the results for r — > 2n -k + ?>it®v t and r — ► 3ir 2tt + it®v t . All errors 
are statistical, except the second errors in the branching fractions, which are systematic. 



Decay mode 




2ir 7r + 37r°z/ T 


3ir 2ir + TT°i' T 


Data (events) 
qq bg (events) 
r bg (events) 
Efficiency (%) 




139.0 ± 11.8 

15.1 ±3.1 

35.2 ±3.4 
1.65 ±0.03 


231.0 ± 18.8 
25.8 ±5.9 
19.4 ± 5.5 

4.45 ± 0.06 


B (HT 4 ) 




2.2 ±0.3 ±0.4 


1.7 ±0.2 ±0.2 


TABLE II. Summary of the results for the intermediate states with an rj meson. All errors are 
statistical, except the second errors in the branching fractions, which are systematic. 


Decay mode 


2tt TT + r]v T 
rj -» 3vr° 


7r 2^r\v T 

7] — > 7r + 7T _ 7r° 


2tt TT + r]v T 
r\ — > 7r + 7r _ 7r 


Data (events) 
qq bg (events) 
r bg (events) 
Efficiency (%) 


32.1 ±6.7 
1.9 ± 1.5 
0.9 ±0.9 

1.28 ±0.05 


15.4 ±5.4 
2 +1 '° 
2.3 ±1.6 
1.48 ± 0.05 


52.4 ± 10.4 
5.2 ±3.2 
2.5 ±2.2 

4.18 ±0.08 


B (HT 4 ) 


2.9 ±0.7 ±0.5 1.5 ±0.6 ±0.3 1.9 ±0.4 ±0.3 



combinations in the u mass region are shown in Fig. |3|. There is also a clear signal at the uj 
mass. In particular, the decay r _ — > 2-k^-k + ujv t is observed for the first time. 

The r) and uj signals are extracted by maximum likelihood fits with a dual Gaussian over a 
polynomial background. The shapes have been constrained to the Monte Carlo expectations. 
The results are summarized in Tables |TJ and 

There are several sources of systematic errors. These include the uncertainty in the 
number of produced r + r~ pairs (1.4%), charged track reconstruction (1% per track), 7r° 
reconstruction (4% per 7r°), efficiency (1-5%) and r migration background (4-12%) estimates 
due to limited Monte Carlo statistics, hadronic background estimates (4-9%), mixtures of 
various intermediate states (1-6%) and branching fractions of rj — > 3tt° and tt + 7t~7t° and 
uj — > 7T + 7r~7r° 0. For the extraction of uj and r] signals, there are also systematic errors 
(5-10%) resulting from the choice of combinatorial background shape and fit region. The 
branching fractions with systematic errors are listed in Tables |-|I|. The results represent 
significant improvement in precision over previous measurements 0. The branching fractions 
for the two decays with uj in the final states are somewhat smaller than the recent calculations 
by Gao and Li 0. 

The results on the six-pion decays can be compared with the isospin symmetry and 
CVC predictions, after correcting for the contributions from the axial-vector current r~ — > 
(37r)~?7z/ r , which also violates isospin conservation with the decays r] — > 3tt° and 7r + 7r~7r°. 
To reduce the uncertainty in the corrections, we use measurements from these decays and 
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TABLE III. Summary of the results for the intermediate states with an u meson. All errors 
are statistical, except the second errors in the branching fractions, which are systematic. 



Decay mode ir 2ti^ujv t 2ir tt + uju t 

Data (events) 54.3 ± 11.4 110.1 ± 18.6 

qq bg (events) l.Otji 3.5±f;§ 

r bg (events) 8.7 ± 3.4 2.3±|| 

Efficiency (%) 1.39 ±0.06 4.06 ±0.14 

B (10" 4 ) 1.5 ±0.4 ±0.3 1.2 ±0.2 ±0.1 



the decay rj — > 77 J/J to obtain the average branching fractions: B(t — > 27r n + r]v T ) = 
(2.4 ± 0.5) x lCr 4 and B(r~ -> tt^ttVt) = (1.5 ± 0.5) x 10~ 4 . This yields the vector 
current branching fractions: By{r~ — ► 27r _ 7r + 37r°z/ T ) = (1.1 ± 0.4) x 10~ 4 and By{j~ — > 
3/T _ 27T + 7r°z/ r ) = (1.1 ± 0.2) x 10~ 4 , corresponding to ~50% and ~65%, respectively, of the 
inclusive six-pion branching fractions. 



The isospin model |LT] classifies n-pion final states into orthogonal isospin partitions and 
determines their contributions to the branching fractions. The partitions are labeled by three 
integers (ni, 712, ris), where is the number of isoscalar subsystems of three pions, n 2 — n% 
is the number of isovector systems of two pions, and ri\ — n 2 is the number of single pions. 
For n = 6 there are four partitions: 510 (47rp), 330 (3p), 411 (3iiu), and 321 (vrpcu), denoted 
according to the lowest mass states. The isospin model imposes constraints on the relative 
branching fractions, which can be tested by comparing the fraction 

B v (t- -> 27r-7r+37rV) 

I 27r-7r+37r° " 



B v (t- -> (6tt)-i/ t ) 

to a similar fraction f3 7r -2 7 r+7r°) where By{r~ —* (67r)~z/ T ) is the sum of the branching fractions 
for the three six-pion vector decays. Figure § shows f 27r - 7r +3 7r o vs. f 37r - 27r + 7r o with our new 
measurement of the branching fractions. The measurement is presented as a line because 
By{r~ — > 7r~57r°z/ r ) has not been measured yet. The result is consistent with the isospin 
expectation since the experimental measurement overlaps with the isospin triangle. Our 
result indicates the 321 (npu) partition is dominant because the decays t~ — > h~2t{ oji> t and 
t~ — ► 2/T~ tt + ujv t saturate the six-pion (vector) decays. 

The CVC hypothesis combined with isospin symmetry predicts that the branching frac- 
tions for both t~ — > 27r~7r + 37r°z/ r and t~ — ► 37r _ 27r + 7r°z/ r are greater than (2.5 ± 0.4) x 
10~ 4 [[L^]. This is somewhat larger than the measured branching fractions for the six-pion 
vector decays. The discrepancy is even more significant if we compare the predicted inclusive 
branching fraction Bcvc{j~ ~ > (67r) _ z/ r ) = (12.3 ±1.9) x 10~ 4 with the sum of the measured 
six-pion vector branching fractions under the assumption that By{j~ — > 7r~57T°z/ r ) is com- 
parable with or smaller than By{j~ — > 27r~7r + 37r°z/ T ) and By{j~ — > 37r~27r + 7r°z/ T ) JT5|| . This 
assumption is consistent with the observation that the six-pion vector decays are saturated 
by intermediate states with an u meson, which implies a small decay width for the 510 (47rp) 
state, the only state that contributes to the decay t~ — > 7r _ 57r°z/ r . The discrepancy might 
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be explained by sizable presence of / = states in the e + e~ annihilation data that should 
be subtracted before calculating the CVC prediction. 

In conclusion, the branching fractions for two six-pion decays have been measured with 
much improved precision. The resonance substructure in the decays has been studied. In 
particular, the decay t~ — > 2it~-K + ujv T has been observed for the first time. Within the 
statistical precision, the decays are saturated by rj and uj intermediate states. The measured 
branching fractions are in good agreement with the isospin expectations but somewhat below 
the CVC predictions. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, the Research Corporation, the Natural Sciences and 
Engineering Research Council of Canada, the A. P. Sloan Foundation, the Swiss National Sci- 
ence Foundation, the Texas Advanced Research Program, and the Alexander von Humboldt 
Stiftung. 
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FIG. 1. Mass spectra of the hadronic systems in the decays (a) r~ — > 2tt~7t + 3tt i' t and (b) 
r~ — > 37r _ 27r + 7r ^ T . The solid histogram is the sum of the signal Monte Carlo and background 
(dashed), which includes the r migration and hadronic (shaded) background. 
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FIG. 2. Three-pion mass spectra of r\ candidates: (a) Af 37r o in r — > 2ir 7r + 37r°z/ T , (b) M n +^- n o 
m r ' -» 2tt-7t + 3tt u t , (c) ikf 7r + 7r - T o m r There are six entries per event in (b) 

and (c). The solid lines shows fits to the data. 
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FIG. 3. Three-pion mass spectra of u candidates (six entries per event) in (a) 
r~ — > 27r _ 7r + 37r z^ T and (b) t~ — > 37r _ 27r + 7r°^ r . The solid lines shows fits to the data. 




FIG. 4. Decay fractions of t~ — > (Qir)~v T . The thick solid line through the origin represents 
the measurement. The shaded area indicates the one standard deviation region, calculated with 
correlated errors taken into account. The triangle bounded by the dots shows the isospin expecta- 
tion. 
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